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ABSTRACT 

o 

(T) ' Differential regularization is applied to a field theory of a non-relativistic charged 

| boson field (ft with X((ft*(ft) 2 self-interaction and coupling to a statistics-changing 



Q_i' U (1) Chern-Simons gauge field. Renormalized configuration-space amplitudes for all 

diagrams contributing to the (ft* (ft* (ft(j) 4-point function, which is the only primitively 
divergent Green's function, are obtained up to 3-loop order. The renormalization 



X 

group equations are explicitly checked, and the scheme dependence of the /3-function 
is investigated. If the renormalization scheme is fixed to agree with a previous 1- 
loop calculation, the 2- and 3-loop contributions to /3(A, e) vanish, and /3(A, e) itself 
vanishes when the "self-dual" condition relating A to the gauge coupling e is imposed. 
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I. INTRODUCTION 



The differential renormalization procedure is a novel method for perturbative calculations in 
quantum field theory in which real space correlation functions are simultaneously regularized and 
renormalized. Neither explicit cutoff nor singular counter terms are required, yet no ultraviolet 
divergences appear. The method has now been applied to many relativistic field theories, and 
a systematic formulation ||, obtained with the use of BPHZ topological analysis of Feynman 
diagrams, has recently been given. 

Field theories with non- relativistic kinematics are also of physical interest, and we report in this 
paper on the application of differential regularization to a theory in two spatial dimensions which 
has been used to describe the anyon Aharonov-Bohm and fractional quantum Hall effects. The 
particular theory [[3] we study has evolved from earlier work It describes a non-relativistic 

boson field with quartic self- interaction and coupling to a 2+1 dimensional abelian gauge field with 
Chern-Simons kinematics. 

It was shown in || that the classical Lagrangian of the theory is scale invariant, but one expects 
this symmetry to be broken by renormalization effects at the quantum level. These scaling properties 
are similar to those of renormalizable massless relativistic theories in 3+1 dimensions, so the theory 
is a natural place to apply a new regularization method. The theory is also chiral so that dimensional 
regularization is problematic. 

The renormalization properties of the theory have been studied at the 1-loop level ||[F| using a 
momentum cutoff to handle ultraviolet divergences. Scale symmetry breaking is found for general 
values of the coupling constants, but scale invariance is restored when the couplings are related by 
the "self-duality" condition which has special significance in the classical theory [[|. 

In this work we use the differential regularization method to renormalize the divergent Green's 
functions of the theory through 3-loop order. This method is based on the observation that bare 
amplitudes which are primitively divergent in momentum space are well defined for separated points 
in real x, t space, although too singular at short distance to possess a Fourier transform. These 
amplitudes are then renormalized by the following two-step procedure: 

(1) Use differential identities (typically containing a new mass scale parameter M) to express 
singular bare amplitudes as derivatives of less singular Fourier transformable functions. 

(2) Define the renormalized Fourier transform and other integrals of the originally singular am- 
plitudes by formal partial integration of derivatives. 

Multi-loop amplitudes are then handled by inserting the differential identities for their primitively 
divergent subgraphs, manipulating differential operators according to step (2) and using new dif- 
ferential identities as needed. The heuristic reason that this procedure is consistent is that the 
ultraviolet singular surface terms which are neglected in rule (2) correspond to singular counter 
terms which occur in traditional regularization methods. The operational test of consistency is that 
the amplitudes obtained obey renormalization group equations (RGE's) in the scale parameter M, 
with /3-functions and anomalous dimension defined through the equations themselves. 

New features occur when these ideas are implemented in a non-relativistic field theory. Because 
of the lack of symmetry between space and time coordinates, differential identities of rather different 
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structure from those for relativistic theories must be found. Due to the combined non-relativistic 
and Chern-Simons kinematics the gauge coupling is not renormalized, so that there is only a (3- 
function for the quartic coupling. Nevertheless this /3-funtion is scheme-dependent beyond one-loop 
order, as expected for a theory with two or more couplings. In differential renormalization, scheme 
dependence appears because of the freedom to introduce different mass scales M, M' in some of 
the differential identities used to regulate the overall divergence of graphs which are not related by 
symmetry. If this occurs in a given order of perturbation theory the /3-function typically depends on 
log M/M' in the next order and beyond. Our calculations show that there are two scale parameters 
which appear at the one-loop level and one more at the three-loop level. The dependence of the 
two- and three-loop /^-function on the ratio of the one-loop scale parameters is obtained. If the 
mass ratio is fixed to make our one- loop result agree with then the two- and three- loop 

/3-function vanishes. The remaining one-loop contribution then vanishes when the two couplings 
satisfy the "self-duality" condition of ||. Thus previous results on the relation of "self-duality" to 
scale invariance in the quantum theory have now been extended through three-loop order. 

Our approach to this problem has been purely field-theoretic; we have simply obtained renor- 
malized real space Green's functions which obey RGE's. We have not compared results with closely 
related quantum mechanical studies of anyons || , although it appears that some difficulties of per- 
turbative quantum mechanical treatments can be resolved using field-theoretic methods. Finally, 
we can only express the hope that our work can be applied to elucidate physical properties of anyon 
systems.^ 

The Lagrangian of the theory under study is 

L = <p* (iD t + ^-D 2 ) + at, l3 (^dtAiAj - A^A^j - A(0*0) 2 (1.1) 
with covariant derivatives 

D t = d t + ieA 

(1.2) 

D = V-ieA 

In the Coulomb gauge the only non-vanishing components of the gauge field propagator take the 
instantaneous form [] 

(Mx)A (0)) =iTi{x) 

- s £ * s (1 . 3) 



1 During the writing phase of this work, references || came to our attention through the bulletin board. 
Three-loop calculations of the finite V, T partition function of the theory at the self-dual point are presented 
in these papers. 

2 We use the notation x = (x, t). 



as required by the physical role of as a statistical gauge field without propagating degrees of 
freedom. The Schrodinger propagator is 

(T0(:r)0*(O)) = iG(x) 

m • > (1-4) 



It satisfies 



id t + G(x) = 5(t)5{x) (1.5) 
\im + G(x) = -i(5(x). (1.6) 



Perturbation calculations may be performed with the propagators ( |1.3|) , (1.4) and the interaction 



vertices from (1.1). The superficial degree of divergence for the different Green functions can 
be calculated by taking into account that the scale dimension of time is twice that of the space 
coordinate. 

Many Feynman graphs vanish due to the special kinematics of the theory. There is no particle 
production so the number of 0-field lines is conserved in intermediate states of any graph. Many 
graphs with internal gauge field lines vanish when one adds the contributions where attachments 
of A and Aj are exchanged. These features of the theory imply that there are no quantum loop 
contributions to the propagators (1.3) and (1.4), the <ft*<j)A vertex and the <j)*<j)AA "Compton" 
amplitude. This means that the only superficially divergent Green's function is the 0*0*00 4-point 
function, and there are no anomalous dimensions in the theory. From the 4-point amplitude one can 
construct 6-, 8-, 10-point functions, etc., using a skeleton expansion without encountering further 
ultraviolet divergences. 

In Section II below we confine our attention to the scalar subtheory of (1.1). We discuss the 
basic differential identity for the 1-loop scalar bubble graph, and its Fourier transform. The convo- 
lution theorem can then be used to compute and sum multiple bubbles and obtain the full off-shell 
scattering amplitude. In Section III, the regularization of the 1-loop graphs of the full theory is 
obtained. Two- and three-loop contributions are considered in Sections IV and V, respectively. The 
verification of the RGE's and the /^-function are discussed in Section VI. 



II. SCALAR SECTOR 

In this section we will study the scalar sector of the theory defined by the Lagrangian (1.1). The 
tree level contribution to the amputated 4-point function is 

r a (a ; j) = —4:i\5(xi — x 2 )<5(xi — x 3 )S(xi — x 4 ) (2.1) 

where we have used the simplified notation 8(x) = S(t)8(x). 

In a non-relativistic theory, the only non-vanishing one-loop diagram is the direct channel bubble 
graph depicted in Fig. lc. The bare amplitude in configuration space is given by 
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T c (xi) = 8\ 2 5(xi — X2)5(xz — X4)G 2 (x3 — x\) 



(2.2) 



where G(x^ — Xi) is the scalar propagator given by eq. (1.4). This amplitude is well-defined for 
separated points, but the factor 1/t 2 is too singular at short distances, yielding a logarithmically 
divergent Fourier transform. 

The bare amplitude can be regulated by the following differential identity 

G 2 req {x) = ~ (idt + F(x) 

(2.3) 



An 2 \ ' Am) 
F(x) = \ogiM 2 t- 



2 0[is\ it/ x 



t 

which is an exact relation away from the origin constructed to incorporate the ideas discussed in 
the introduction. 

i) The differential operator acts on the function F(x), which is less singular than G 2 (x) and has 
a well-defined Fourier transform. Using Gaussian integration it is straightforward to obtain 

~ 7T log(k 2 /4m -u- xe^/M 2 } 

l<(U,k.) = — — ; , (2.4J 

m k^/Am — uj — le 

where 7 = 1.781 .. . is Euler's constant. 

ii) The coincident point singularity of (2.2) is defined as in distribution theory using formal 
partial integration of derivatives in any integral. For the Fourier transform, this prescription gives 
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f 2im\ 
T c (u, k) = log 

7T 



Am j M 2 



(2.5) 



iii) Most important is the fact that the parameter M, which is required for dimensional reasons 
in ( p.3|) , plays the role of the renormalization group scale parameter. To see this we can compute 
the effect of a variation of M 



, r 9 „o / >. im 2 



id t + 



Am) t 



(2.6) 



2n 



where we have used ( |1.6| ) and the fact that the logiM 2 £ factor in (2.3) multiplies a Schrodinger 



propagator with mass 2m. Comparing (2.1, 2.2 , 2.6), we see that the effect of a change in M can be 
compensated by a change in the coupling constant. Indeed from the renormalization group equation 



M ii T 



-/3(A) 



d_ 



X i 



(2.7) 



with taken as the sum of the tree and renormalized bubble amplitudes one finds the one-loop 
/3-function 
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/3(A) = — (2.8) 

7T 

iv) The factor of i in \ogiM 2 t is required to obtain a real analytic amplitude in momentum 
space with non-zero imaginary part for uj > k 2 /4m as required by unitarity. Alternatively, as in 
the relativistic case [1], one could start in the imaginary time formalism, t = —it, in which all bare 
amplitudes and the differential identities used to regulate them are real functions of x and r. Then 
logiM 2 t is automatically obtained for real Schrodinger time. 

In the scalar theory it is easy to sum multi-bubble graphs and obtain the full scattering ampli- 
tude. These graphs are convolutions in real-space and the renormalized amplitudes in p-space can 
be defined as products of the renormalized 1-bubble amplitude (2.5). The scattering amplitude is 
then 



CO 



A(u;,k) = -4a£(^f>,k)J . (2.9) 
In the cm. frame, k = 0, one obtains 

A(u, = r -, rr-. 2.10 

"2T g M 2 

The above remarks are intended to show how the differential renormalization procedure repro- 



duces previous results on the scalar sector [10,11], obtained in most cases with a momentum cutoff 



and renormalization at a momentum scale p 2 = fi 2 (here p = Pi = — P2? where pi and P2 are the 
momenta of the incident particles). Since the on-shell cm. total energy is given by uo = p 2 /m, 
our renormalized amplitude can be brought into full agreement with previous work if the scale 
parameters are related by /1 2 = mM 2 /^. 



III. ONE LOOP ORDER 

In this section we shall illustrate how differential regularization works for the full theory to 
one-loop order. The tree graph of Fig. lb does not require regularization, but for completeness we 
give its x-space amplitude, 

e 2 d 
T h (xi) = — 5{xi - x z )Ti(xi - x 2 )tt-S(x 4 - x 2 ) (3.1) 
m obq 

The one-loop graphs are shown in Fig. lc-lf. In addition to the diagrams explicitly drawn 
one must add, as appropriate in each case, i)the exchange of A and Aj vertices, ii)the opposite 
orientation of each triangle with Aq, Aq, A 2 vertices, iii)the time-reflected graph with initial and final 
states interchanged, iv)exchange of x 3 and x 4 as required by Bose symmetry. It would be tedious 
to discuss these permutations explicitly, but it is important to include them when the RGE's are 
discussed in section VI. 
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Power counting leads us to expect that all one-loop graphs for the 4-point function are ultraviolet 
divergent. However, as we shall show, the only 1-loop graph which actually requires renormalization 
(apart from the bubble already discussed in section II) is the seagull graph shown in Fig. Id. We 
therefore begin with this graph whose bare amplitude is given by the following expression, 



e 4 

T d (xi) = —5(x 2 - x A )G(x 3 - xi)Ti(x 2 - x 1 )T i (x 2 - x 3 ). (3.2) 
m 

Note that in this expression there is a factor 

0(*3-ti)<*(*3-*i) (3-3) 

where the 9 and 5 function come from the boson and gauge field propagator respectively. We then 
use that, 

9(t)5(t) = l -5(t) (3.4) 

a property that can be deduced from the convolution theorem. Using this together with the value 
at t = of the boson propagator (|1.6| ), we obtain, 

ie 4 1 

rd (^) = - o ^ 5 ( x * - xt)8(xi - x 3 )8(t 3 - £2)7 (3.5) 

8mir z a 2 (x 2 — xi)2 

Notice that as the boson propagator at t = is proportional to a delta function, the triangle 
"collapses" to an effective gauge field bubble. The function 7^3^32 is too singular at short distances 
and as a result, its Fourier transform is logarithmically divergent. In the spirit of differential 
regularization, we use the following identity, valid at all points except x = 0, 

-I = I V 2 log 2 M 2 x 2 (3.6) 
x^ 8 

Note that we have introduced a different mass scale from that in the identity ( |2.3| ) associated 
with the bubble graph. A priori, one can use independent scales for different graphs which are not 
related by a Ward identity. Different choices of M's correspond to different renormalization schemes 
and in section VI we shall see how this affects the /3-function. 

After using equation fl3~6[ ), we obtain the regularized form 



• 4 

le 



T d reg (x t ) = - -^—^ 5(x 2 - x A )6( Xl - x 3 )5(t 3 - t 2 )W{ log 2 M s 2 ( Xl - x 2 ) 2 (3.7) 

which has now a well defined Fourier transform, since it can be shown using the same technique as 
in the Appendix of Q that 

d 2 xe^ log 2 M 2 x 2 = -% log ^ (3.8) 

From this expression we can calculate the contribution of this graph to the 4-point function in 
momentum space. After taking into account the graphs obtained by permutation, we have in the 
center of mass frame, 
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r D (p) = -^iogf^^) (3.9) 

2ira z | singly 



which coincides with the result in ref. [|7| if the identification [i = M s = is made. 

Comparing (|3.9| ) and the momentum space amplitude of the bubble graph fl2.5|) in the cm. frame 
we also see that the renormalization scheme adopted in Q corresponds to the following relation 
between scales 

M 2 = Pl M 2 = ^M 2 (3.10) 

We end this section by showing that the remaining one-loop graphs do not require renormaliza- 
tion. It can be easily shown that contribution of the triangle graph shown in Fig. le cancels against 
the one coming from the graph with the < A Ai > propagator running in the opposite direction. 
Of course, any higher-loop diagram that contains this one as subdiagram will automatically vanish, 
therefore they are not included in Fig. 1 and we shall omit them from our discussion. 

We shall show next that contributions of the box graphs are finite. Let us consider as an example 
the diagram shown in Fig. If, which is given by 

T f ( Xl ) = -^G(x 3 - x x ) ( / d G {x A - x 2 )) %{x 2 - x 1 )T j (x 4 - x 3 ) (3.11) 
m z \Cx 2 c>X4 / 

Because the gauge field propagator is transverse, 

V x -T(x-y) = (3.12) 
one can immediately extract the derivatives as total derivatives and rewrite the amplitude as, 

T f (xi) = — — r- — — r [G(xs - xi)G(x A - x 2 )Ti(x 2 - x 1 )T j (x i - x 3 )} (3.13) 

The Fourier transform of this amplitude can be expressed now as a finite loop momentum 
integral and the total derivatives are simply factors of external momenta. The other box graphs 
with A <-» Ai permutations can be treated similarly. 



IV. TWO LOOP ORDER 

At two-loop order there are two divergent graphs: the double bubble shown in Fig. lg and the 
"ice cream" graph shown in Fig. lh. 

A regularized expression for the double bubble is obtained immediately using the regulated 
expression (|2.3| ) for each bubble subgraph, 

r 9 reg ( Xi ) = 16iA 3 5(xi - x 2 )6(x 3 - x 4 ) J d 3 uG 2 reg (x 3 - u)G 2 reg (u - x x ) (4.1) 

This integral can be explicitly computed and we obtain 
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7T 



( v 2 \ 

5(x 1 - x 2 )5(x 3 - x 4 ) I id t3 + — - I S(x 3 - x x ) 



(4.2) 



where 



0(t) jmx 

5(a;) = e * 



log; 2 iM 2 t - 



7T 



(4.3) 



Now, let us turn to the regularization of the ice cream graph. The bare amplitude is given by 



T h (x t 



e 4 A 



2nrn 2 a 2 



8(xi — X2)S(ts — t^)G(x3 — X2)G(X4 — X\)- 



x 3 - x 4 ) 2 



(4.4) 



For simplicity, let us set X\ = and introduce coordinates x, y and t such that X3 = (x, t) and 
24 = (y, i). The singular part of the graph is then contained in 



/ /l (x,y,t) = G(x,t)G(y,t)- 



x-y 



(4.5) 



We begin by regularizing the divergent seagull subgraph (see Eq. ( |3.6|) ) 

/ fe (x, y, t) - / fc (x, y, t) = ic(x, t)G(y, t)(V, - V,) 2 log 2 M 2 (x - y) 2 
After manipulating derivatives, this expression can be rewritten as 



(4.6) 



Ax,y,i) 



1 

32 



-(v> 



log 2 (M 2 (x-y) 2 )(Vx - V y )G(x,t)G(y,t) 



+ log 2 M 2 (x-y) 2 (V a; - V y yG(x,t)G(y,t) 



(4.7) 



The first term in this expression, being a total divergence, is finite by power counting. In order 
to improve the second term, note that the product of two non-relativistic propagators satisfies the 
identity 



(Vx- V y ) 2 G(x,t)G(y,t) 



(V* + V y f + &mdt 



G(x, t)G(y, t) - 4zm5(x)5(y)5(t) (4.8) 



To regularize ( f4.7| ) we drop the last term in this identity, because it simply gives a singular surface 
term in the Fourier transform of f h . In this way, we obtain the regularized expression 



32 



log 2 M 2 (x-y) 2 (Vx - V y )G(x,t)G(y,t) 



+ 



(Vx + V y ) 2 + 4imd t 



log 2 M 2 (x-y) 2 G(x,t)G(y,t) 



(4.9) 



This can be inserted in ( |4.4| ), restoring the original space-time points, to obtain a regularized 
form for the amplitude T h (xi). 

We shall show next that the rest of the two-loop graphs in Fig. 1 are indeed finite. The bare 
amplitude for diagram i is 



8 



Fix,, 



le 



2An 2 a 2 



_d_ 

<9x? 



■G(X4 - x 2 ) 



G(x 3 - xi)Ti(x 2 - xi 



x 4 - x 3 ) 



(4.10) 



We first regularize the subdivergence x 3 — > x 4 by using the identity 



ie 



24n 2 a 2 



_d_ 

dx.1 



G(X4 - x 2 ) 



G{x 3 - xi)T<(x 2 - x ± )-V 2 log 2 M 2 (x 4 - xg) 1 



(4.11) 



Using now the property fl3.12Q , the overall divergence can be regularized by writing the diagram 
as total derivative, 



Fix*) 



ie 







24vr 2 a 2 dxi 



G{x A - x 2 )G{x 3 - x 1 )T l (x 2 - x^-V 2 log 2 M 2 (x 4 - x 3 ) ; 



(4.12) 



Diagram i is now regularized. Nevertheless, let us point the following fact. Though this diagram 
depends on M s , the sum of this diagram plus the one where the < AqA, > line runs in the opposite 
direction does not. This can be understood in the following way: as easily seen from ( 4.12| ), the 
scale derivative of F l is proportional to T e , and we have already shown that this diagram is zero 
once the corresponding Aq <->■ Ai exchange is taken into account. We then see that although a scale 
dependence is introduced in intermediate steps to properly handle sub divergences, this dependence 
can disappear in the final amplitude. As we shall see, a similar situation occurs in several three-loop 
diagrams. 

The last two- loop diagram, Fig. lj, can be shown to be finite. After some manipulation of 
derivatives, the amplitude can be written as 



THX; 



d d 



m 3 <9xn dxh 



Ti(x 2 - xi)T k (x A - x 3 ) 



d 
dxi 



G(x 3 — v)G(x4 — u)G{u — x 2 )G{v — x\)Tj(u — v) 

(4.13) 



The Fourier transform of this amplitude is a two-loop energy-momentum integral which is convergent 
by power counting because it contains explicit factors of external momenta. 



V. THREE LOOP ORDER 

There are eight diagrams contributing to the 4-point function at three loops (see Fig. 1) although 
we shall show that only four of them are divergent, namely diagrams k-n. To regularize three- loop 
bare amplitudes we first identify the divergent subgraphs and treat them in the same way as we 
have done at lower orders. The regularization of overall divergences (when all the external points 
come close) may require new DR identities. Our main concern is to compute the contribution of 
every diagram to the RGE's (see section VI) and thus in some cases we won't give the explicit 
closed form of the diagram, since it is not needed, but only an integral representation that has an 
ultraviolet finite Fourier transform and from which we can obtain the relevant contributions of the 
diagram to the Callan-Symanzik equations. 
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Let us start with the triple bubble diagram. We regulate each bubble subgraph using eq. (|2 
and we obtain 

^reg( x i) = -32A 4 <5(xi - x 2 )5(x 3 - x 4 ) / G 2 eg (u - xx)G 2 reg {v - u)G 2 eg (x 3 - v) (5.1) 

Ju,v 

This expression is the convolution of three factors which have a finite Fourier transform, therefore 
so does the full result, and we shall find the contribution to the RGE's explicitly from fl5.1| ). 

Diagram 1 has also the structure of a convolution of two subdiagrams we have already regulated, 
namely the one bubble graph (|2.3| ) and the "ice cream cone" ( f4.9|) . Using these we obtain 



reg 



iX 2 e 



2^4 



m(ira 



;S(Xi - X 2 )S(t 3 - U) 



X 



J dtd 2 uG 2 reg (u - x x \ M)tf eg (x 3 - u, x 4 - u, t 3 - t; M s 



(5.2) 



where, as we have explicitly indicated, the mass scales used in regulating the sub divergences are 
different. Again, ( |5.2j ) is an adequate regulated expression, since each factor in the convolution has 
a well defined Fourier transform, and thus it is all we need to find the contribution of the diagram 
to the Callan-Symanzik equations. 

The regularization of diagram m is more involved. The bare amplitude is 



T m (x,) 



64m 2 



(5.3) 



where 



f m (x t 



(X 2 - Xi) 2 (x 4 - X 3 ) 2 



G(x 3 — xi)G(x4 — x 2 ) 



(5.4) 



After relabelling points Xi — > x, x 2 — » y, x 3 — >• z, x 4 — >• 0, t\ = t 2 — >• 0, t 3 = i 4 — >• t and regulating 
the two divergent seagull subgraphs by using Q3.6|), we obtain the partially regulated expression: 



/'" (x. y. z. / ) = ^ 2 G(z - x, t)G(-y, t)(V x - V y ) 2 log 2 M s 2 (x - y) 2 V 2 log 2 M 2 z 2 



(5.5) 



Now we integrate by parts the differential operator (V x — V y ) and we use the relation (|4.8| ) satisfied 
by the product of two scalar propagators. Note that when we insert the exact identity ( f4.8p in ( |5.5| ) 
we have to keep the term that contains the ^-functions because it is not a purely local contribution 
(as it was in the "ice cream cone" amplitude (|4.8|) , where we dropped it), but it multiplies a factor 
which requires further regularization. 
We then rewrite (5.5) as 



/ m (x,y,z,t) 



16 2 



(Vx - V 5 



log 2 M 2 (x - y) 2 (Vx - V y )G(z - x, t)G(-y, t) V 2 log 2 M 2 z 2 



+ 



(Vx + V y ) 2 + 4imd t 



2„2 



log 2 M 2 (x - y) 2 G(z - x, t)G(-y, t) V 2 log 2 M 2 z 



\og z M 2 x 2 



+32im5(t)8(z - x)5(y) 



(5.6) 
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Note that once we have regulated the sub divergence z ~ by using ( |3.6| ), we do not need to 
integrate by parts the operator V 2 as we did with (V x — V y ) 2 . The reason is that the first two 
terms are already total derivatives with respect to external points and, since the diagram is only 
logarithmically divergent, these terms have a finite Fourier transform, defined by partial integration 
neglecting surface terms. The remaining term is regulated using the new DR identity 



log z M 2 x 2 



^V 2 log 4 M 2 x 2 



(5.7) 



Thus the final regulated amplitude is 

,4 \ 2 1 



r m (x- 



2m 



lira 



-5{t 2 - t x )S{U - h) 



x (V 



log 2 M 2 ( Xl - x 2 ) 2 (Vi - V 2 )G(x 3 - x l )G{x A - x 2 )V 2 log 2 M 2 (x 3 - x 4 ) 



+ 



(Vi + V2) 2 + 4imd t3 
2im 



log 2 M 2 ( Xl - x 2 ) 2 G(x 3 - xi)G(x 4 - x 2 )V 2 log 2 M 2 (x 3 - x 4 ) 



5(t 3 - tJSfa - x 3 )5(x 2 - x 4 ))V 2 log 4 M 2 ( Xl - x 2 ) 2 | 



The bare amplitude of diagram n is given by 

!-<*) = ^ F 2 
m \7ra 

where 



S(Xi — X 2 )S(x 3 — Xi)l(x 3 — X\ 



I(x 3 - xi) = [ G(u - xi)G{v - x l )G(x 3 - u)G(x 3 - v)8(t u - t v )- — — 

Ju,v (U — Vj' 



(5. 



(5.9) 



(5.10) 



This is the first diagram which involves an internal integration which is not a convolution. Indeed 
this integral contains a logarithmic subdivergence for u ~ v. Our strategy is to regulate the 
singular factor l/(u — v) 2 using ( |3.6|) and then to perform the integral over internal points using 
formal partial integration. After this is done we will cure the overall divergence when x 3 ~ x\. 

We relabel points as x\ — * 0, x 3 — > (x, t), t u — t v — > t' . After regulating the divergent seagull 
subgraph by using (|3.6| ) we integrate by parts the laplacian operators. Total derivatives with respect 
to internal points can be dropped and we obtain the partially regulated form 



J(x, t) = — [ dt'd 2 ud\ log 2 M s 2 (u - \ 
32 J 

(V u -Vv?G(u,t')G(v,t') 



(5.11) 



G(x-u,t-OG(x- v,t-t') 



+ G(u, f)G{v,f) 



(Vu - Vv) 2 G(x - u, t - t')G{x -v,t- 1') 



+2 



(Vu- Vv)G(u,t / )G(v,t / ) 



(Vu - Vv)G(x - u, t - t')G(x - v, t - t') 



The first term, which we shall denote ii(x, t), can be easily regulated by means of the DR identity 
(fOI). The delta term can be dropped because its effect is cancelled by the counterterm which 
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cancels the corresponding singular surface term in the ice cream cone subgraph. After partially 
integrating the differential operator and replacing derivatives with respect to internal points by 
derivatives with respect to x, t, we can rewrite I\ (x, t) as 



/i(x,t) = - — (Aimdt + V 



(5.12) 



x j dt'd 2 ud 2 v\og 2 M 2 (u- v) 2 G(u,t')G(v,t')G(x- u,t - t')G(x- v,t - t') 



The integral in ( 5.12 ) is now finite by power counting and it can be done using Gaussian integration, 
leading to the regulated result: 



h{x,t) 



m 



128tt : 



:(Aimd t + Vl 



6(t) i. 



x 



log 2 iM 2 t - 2(logpi + 2) \ogiM 2 t + k] } 



(5.13) 



with k = log 2 pi + 41ogpi + 8 — n 2 /6 and p 1 = mj/A. 



Note that ( 5.13 ) has a finite Fourier transform, since ii(x, t) was logarithmically divergent and it 
contains total derivatives with respect to external points. It can be shown that the second integral 



in (|5.11) is equal to ii(x, t), so we now focus on the last and most involved one 



J 3 (x, t) = ~J dt'd 2 nd\ ^-^- log 2 M 2 (u - v) 2 G(u, t')G(v, t')G(x -u,t- t')G(x - v, t - t') 

(5.14) 



It is convenient to rewrite 



1 



1 /l 1 

+ 



(5.15) 

t'it-t') t \t' t-t'J K J 

so J 3 (x, t) becomes the sum of two terms, which can be shown to be equal using a change of the 
integration variables. Then we use the identity 



m 2 (u - v) 2 



6(t - 1') 



(t - 1') 3 



g 2(t-t') 



^[(u-xJ^Cv-x) 2 ] 



{Aimd t + V 2 )fc^e^ [(u - x)2+(v - x)2] 



t-f 



(5.16) 



to obtain 



/ 3 (x,t) 



£) 4 1(4^ + V 2 )^) 



x / dt' [ rf 2 urf 2 vlog 2 M 2 (u-v) 2 



(5.17) 



io J t' 2 t-t' 

We can easily perform the integral over internal points and we find the partially regulated form 



J 3 (x,t) 



m 



32n 2 t 



(Aimdt + V 2 ) <6(t)e 



- log 2 tM 2 t - (logpi + 2) logzM 2 * 
f- ^log 2 pi + 21ogpi + 4- 1X — 



(5.18) 
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which still contains an overall divergence when x ~ 0, t ~ 0. To regulate it, we integrate by parts 
the time derivative so that we are left with a term which is a total external derivative and therefore 
has a finite Fourier transform, plus three remaining singular terms that can be easily regulated by 
using (|2.3| ) and two new DR identities, namely 



^-e^log^M 2 * 



6{t) 



I- ' ° n — 1 \~ " 'Am) t 

The final regulated form of 7 3 (x, t) is then 



idt + i^) l^ e -T^ i og »+i %M l t 



m 



: {4imd t + V 



6{t) 



-e * 



x 



32~-' ' * ' x { t 

I log 3 iM 2 s t - -(log Pl + 1) log 2 *M 2 t 
o 2 



+ -log Pi +!ogPi + 2 



7T" 

12 



logiM 2 t 



71= 1,2 



(5.19) 



(5.20) 



Note the presence of a new mass scale M c in (|5.20|) , which appears when regularizing the last term 
in (|5.18|) . It is easily seen that the choice of this mass parameter does not modify the /^-function at 
the three-loop level, so we choose its value such that it cancels the local term (last one) in ( |5.13|) 
and it simplifies the expression of the renormalized amplitude. Restoring the original space-time 
points we obtain 



Y n (r- 



i I mXe 
8 ' 
1 



— J 5( Xl - x 2 )5(x 3 - x 4 ) (id t3 + ^ 



#(t 3 — ti) ' m < x 3- x i) 2 

— -e *3-<i 

U - U 



log 3 iM 2 (t 3 - h) - log Pl log 2 iM 2 (t 3 - k 



+ (log 2 Pl - yjlog«M 2 (t 3 -ti) 



(5.21) 



We will show now that the remaining diagrams of Fig. 1 although superficially log divergent are 
instead finite. The contributions of diagrams o, p, q are given by, 



e 4 \ 2 d 



{Tjfa - z 3 )*(*i - t 2 )V 2 log 2 M 2 ( Xl - x 2 ) ! 
/ G(u — xi)G(v — x 2 )G(x 3 — u)G(x4 — v)Ti(y — u) > 

J u,v \ 



32 \ rrrna I <9xi 







e 4 \ 2 d d 



" "" ' 32 \mira J i)x'., i) x \ 



[Ti(x 2 - x 1 )T j (x 4: - x 3 )(Vi + V 3 



x J dtd 2 ud 2 v log 2 M 2 (u - \) 2 G(u - x\)G{v - x 2 )G(x 3 - u)G(x A - v)\ 



A 



reg 



Xj ) 



) 8{xx - x 2 )^j {Tj{x± - x 3 )(Vi + V 3 



16 \mixa 



(5.22) 



(5.23) 



x J dtd 2 ud 2 v log 2 M 2 {u- vfGiu-x^Giv-x^Gixz-^Gix^-v)} (5.24) 
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As in previous cases, the sub divergences have been regularized and derivatives manipulated in order 
to make the integral over internal variables as well as the Fourier transform respect to external points 
well defined. As in the case of diagram i, although each of these diagrams depends on M s , this 
dependence disappears when we take into account the graphs obtained by the permutation A <-> Aj. 
It is easy to prove that the scale derivative of each of these graphs is proportional to a graph that 
contains graph e as a subgraph and as a consequence the scale derivative of these graphs is zero. 

We conclude with the triple box (diagram r). The contribution of this diagram to the amputated 
4-point function is given by 

d d 

r ( x i) = -[ — ) a „i o_j T *( x 2 - Xi)Tj(x 4 - x 3 ) 




x / G(u - xi)G(y - u)G(x s - y)T k (v - u)Ti(z - y) 



'u,v,y,z 



G(x, -z)-^ G{z -v)-^ G{v - x 2 ) ) (5.25) 

Because of the "transversality" of the gauge field propagator, all derivatives inside the integral can 
be rewritten as derivatives with respect to external points and brought to the left to obtain the 
regulated form 

/e 2 \ 4 d d [ d d 

x / G(u-x 1 )G(y-u)G(x 3 -y)T k (v-u) 

Ju,v,y,z 

- y)G(x A - z)G(z - v)G(v - x 2 )} (5.26) 

Away from coincident external points, one can see that the integrals over internal points are now 
finite by power counting. Furthermore, the amplitude ( |5.26j ) has a well defined Fourier transform 
since it contains two total derivatives with respect to external points. 



VI. RENORMALIZATION GROUP EQUATIONS 

In this section we shall show that the renormalized 4-point function satisfies the RGE equation 
and we compute the /3-function to three loops. We recover the one loop result of [[/j], namely the 
zero of the /^-function for the self-dual case, and we find that the relation between the couplings 
for which the /3-function vanishes is scheme dependent beyond one loop. Scheme dependence arises 
from the two scale parameters M and M s introduced initially at one-loop order, and renormalized 
amplitudes depend on the parameter p = M^/M 2 . As explained below fl5.20| ), an additional scale 
parameter in the regularization of graph n was fixed to simplify its renormalized amplitude. This 
scale parameter affects the RGE's beginning at four loops. 

The RGE equation for the 4-point correlation function is 
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Treg = (6.1) 



Recall the anomalous dimension term is absent because there is no field renormalization in this 
theory. In order to check that (|6.1|) is satisfied, we explicitly compute the scale derivatives of 
the renormalized amplitudes obtained in the previous sections and compare order by order in 
perturbation theory with j^T reg . The contribution to the /3-function in each order can be uniquely 
identified from the purely local term, in which acts on the tree contribution (f2.1p. Non-local 
contributions then provide an independent consistency check, since they must appear with the 
correct coefficients to test that sub divergences have been correctly handled by our method. 
The /9-function can be written as a series expansion in the couplings A, e 

e) = M 2 + he 4 + a 2 \ 3 + c 2 Ae 4 + . . . , (6.2) 

where the subscript in the coefficients indicates the loop order. In a theory with more than one 
coupling constant the /3-function is renormalization scheme independent through one loop order but 
not beyond. We will use capital letter superscripts to denote that the amplitudes include now not 
only the contribution of the particular diagram drawn in Fig. 1 but also the contributions of graphs 
related by Bose symmetry as well as time reflection. All permutations of A$, Ai and Af vertices are 
also added. 

We have already computed the one loop /3-function in the purely scalar subtheory (see 
while from ( p.7|) one can immediately obtain the scale derivative of the seagull graph 

d is 
M dM T ° g = " " X2)5{ - Xl ~ Xs) ^ Xl ~ Xa) 

Y A (6.3) 



4:iima 2 \ 
Thus 



ai = - (6.4) 



7T 



h = — - 2 (6.5) 



M^-(rZ q + rZ (l ) = o (6.6) 



and we see that 

d 

when 

4m 2 a 2 ^ 

showing that scale invariance is recovered for the self-dual point 0. 

At two- loop order the /3-function is renormalization scheme dependent, as we shall see. The 
scale derivative of the double bubble ( flip is trivially found to be 
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8 

M——T G 

dM re9 



2mA, 



7T 



-.c 

reg 



(6.8) 



The case of the "ice cream cone" diagram is considerably more complicated. Although it is true 
that the regularized amplitude has a well defined scale derivative, we must be able to recognize in 
the result the amplitude of graphs that have already appeared at one loop order. However, if we 
calculate the scale derivative of f^ eg from fl4.9| ) we find 



\I—f h 

~ - t J re 



dM' 



reg 



;x,y,t) = --(V 
1 



logM s 2 (x- y ) 2 (Vx - V y )G(x,t)G(y,t) 



(V> 



+ 4imd t 



G(x,t)G(v,t)logM 2 (x-y) 



(6.9) 



which is difficult to interpret immediately. One can show that ( |6.9|) and the regulated bubble (^ 
differ at most by a finite local term independent of M. In order to determine this possible /3-function 
contribution, we compute the Fourier transform of ( |6.9| ) (see appendix) and compare it with the 
regularized expression for the bubble. The result is 



M 



_d_ 

dM 



/* (x,y,t;M.) = 27r5(x-y) 



G 2 



reg 



; X ,t;M) + ^log^(t)5(x) 
4tt pi 



(6.10) 



where p = M 2 /M 2 and p\ = m/y/4. Note that we used M s to regulate the subdivergence of 
the graph, but the scale derivative involves the one loop bubble. The RGE requires that the 
renormalized bubble amplitude appear with the assigned scale parameter, i.e., M. 

Combining (|6.10|) with the "ice cream cone" amplitude given by Eqs. (|4.4j ), fl4.5| ), and adding 
the contribution of the diagrams related by Bose symmetry and time inversion, we obtain 



M 



d 

—T H (M 

dM re ° [ s 



lTxm\a 2 rea 



(2vra) 



log— r 

p 



(6.11) 



It is easy to check that 
loop contributions: 



and (|6.11|) are consistent with the one loop result and yield the two 



a 2 = 



c 2 



Pi 



(6.12) 
(6.13) 



One can see explicitly that the coefficient c 2 depends on the mass ratio we have chosen, p = M 2 /M 2 . 
This means that the point for which the /3-function vanishes is also scheme dependent, and in general 
there will be higher order corrections to the one- loop relation (|6.7| ). However, in the renormalization 
scheme defined by 



M 



pl M 2 = ^M 2 
4 



(6.14) 



we have p — p\ and thus there is no contribution to the two-loop /^-function. As a consequence, 
in this scheme the relation ( |6.7| ) between the couplings which gives a finite theory is valid through 
two loops. 



16 



The scale derivatives of the three loop amplitudes are found to be 



m Jm T "° = "^ r - (fU5) 

"m^"- M) = " v W + d^ r - (M) + pi log p rS » (M) (6 ' 16) 

r> 4 

M ^7 V fea = -7-^^ T rea (6-17) 

where we have specified the mass parameter only when the mass scale used in the regularization 
of the three loop diagram is not the same as the mass scale of the graph which appears in its 
scale derivative. In general, it is immediate to obtain these results from the regulated amplitudes 
discussed in section V and only the scale derivative of diagram m requires some further study. As 
in the case of the "ice cream cone" (diagram h), one can show by direct computation that the 
difference between Mttj-tT™ (see ( |5.8|) ) and the regulated amplitude of diagram h is purely local, 



dM reg 

finite and independent of M. Fortunately, to determine this possible local term we only need to 
compute the Fourier transform of a semi-local piece that does not appear in the regulated form of 
the "ice cream cone" , and we easily find that there is no such local term. 

When we substitute the scale derivatives in the RGE (|6.1|), we find the three-loop contributions 
to the /3-function 

a 3 = (6.19) 
b 3 = (6.20) 

47r{2ira) z pi 

Note that if we choose Ml = p\M 2 the three loop contribution to the /3-function vanishes and thus 
the theory is finite at the self-dual point defined by (|6.7| ). 
To conclude we summarize the result for the /3-function 

/3(A, e) = —A 2 - -^e 4 + * log ^Ae 4 - log 2 ^A V (6.22) 
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APPENDIX: 



In this appendix we shall calculate the scale derivative of the ice cream graph h. In order to do 
so, we shall calculate its Fourier transform and show how it is related to the bubble graph c. Let 
us start by recalling that 



R(^y,t) = M^f r h eg 



1 



(v> 



logM s 2 (x - y) 2 (Vx - V y )G(x, t)G(y, t) 



(V x + V y ) 2 + &mdt 



G(x,t)G(y,t)logM 2 (x-y) 



(Al) 



It is convenient to define variables u = (x + y)/2 and v = (x — y)/2. After manipulations of 
derivatives and use of (|1.4j ), the Fourier transform of R can be shown to be 

i?( Pl , p 2 , u) = 4{*p 2 • f 3 + [-p 2 - p 2 + 4mu]f 2 } (A2) 

where pi = p^ + p y , p 2 = p x - P y , and 

f 2 = - [ d 2 ud 2 vdtlog4M 2 v 2 G'(v / 2u,t)G(v / 2v,t)e- ipi U - ip2 - v+ ^ 
8 J 

f 3 = _ f d 2 ud 2 vdt(V„log4M 2 v 2 )G'(v / 2u,t)G'(v / 2v,t)e- ipi U ~ ip2 - v+ ^ 



(A3) 
(A4) 



The calculation of f 3 is straightforward. A conventional non-relativistic loop integral gives the 
result: 



. m , / —Amu + p? 

P2 • f 3 = J ^g — ; — 2 

4 V — Amuj + Pi + p 2 



(A5) 



The loop integral representation of / 2 is infrared divergent because it contains the factor 1/g 2 
which is the Fourier transform of the logarithmic term. Instead we calculate / 2 by the following 
3- step technique: 

i) a trivial calculation of the transform of G(v2u, t) with respect to u and t, 

ii) a not-so-trivial calculation of the transform of log4M 2 v 2 G(\/2v, t) with respect to v and t, 

iii) a final convolution integral over oj to obtain the net result. 

Step ii) is done by considering the function (v 2 )°G(a/2v, t), calculating its Fourier transform 
and taking the derivative with respect to a at a = to obtain the logarithmic factor. Several 
special function formulae |T2|] are required for this task. The final result is 



h 



im 



8[-Pi - Pi + 4mcj] 



log 



I" 



-Amu + p 2 



r 



(-Amu + p 2 + p 2 



2\2 



(A6) 



Then, 



mi 

R(Pl,P2,t) = — lo§ 



-UJ 



1 \ 2 

p| \ Twy 



Am / AMI 



(A7) 
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Now, using (|2.2| ) and comparing with the regulated amplitude of the bubble in momentum space 
( pT5l) we get the relation 

R(p x ,p y ,t) = 2nG 2 (p x + p y ,uj;M^j (A8) 

where pi = 777,7/4 and p = M 2 /M 2 , which is the Fourier transform of (|6.10|) . 

One should note that the regularized amplitude of graph h, see (|4.9[), differs from its scale 
derivative only by the exponent of the logarithmic factor. The loop integral representation of its 
Fourier transform is also infrared divergent. However the transform is well defined, and a technique 
similar to that described above is required to calculate it. 
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FIG. 1. Figure 1. Diagrams which contribute to the 4-point function T(xi, X2, £3, £4) through three-loop 
order. Dots represent the A{ vertex with derivative coupling. Diagrams related by permutations must be 
added as explained in section III. 
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